Wallis-Ramanuj an-Schur-Feynman 

T. Amdeberhan, O. R. Espinosa, V. H. Moll and A. Straub 



Abstract 

One of the earliest examples of analytic representations for tt is given 
by an infinite product provided by Wallis in 1655. The modern literature 
often presents this evaluation based on the integral formula 



2 r°° dx _ 1 /2n 



In trying to understand the behavior of this integral when the integrand 
is replaced by the inverse of a product of distinct quadratic factors, the 
authors encounter relations to some formulas of Ramanujan, expressions 
involving Schur functions, and Matsubara sums that have appeared in the 
context of Feynman diagrams. 



1 Wallis's infinite product for tt 

Among the earliest analytic expressions for tt one finds two infinite products: 
the first one given by Vieta [H] in 1593, 




and the second by Wallis [22] in 1655, 



2_l-3 3-5 5-7 7-9 

^ ~ 2^ ' 4^ ' 6~6 ' 8~8 "' ' 

In this journal, T. Osier [15j has presented the remarkable formula 



(Ll) 



2\ 




2P+^n - 1 2P+^n + 1 



n=l 



2P+h 



2P+h 



where the nth term in the first product has n radical signs. This equation 
becomes Wallis's product when p — and Vieta's formula as p — > oo. It is 
surprising that such a connection between the two products was not discovered 
earlier. 
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The collection [T] contains both original papers of Vieta and Wallis as well 
as other fundamental papers in the history of tt. Indeed, there are many good 
historical sources on tt. The text by P. Eymard and J. P. Lafon [5] is an excellent 
place to start. 

Wallis's formula is equivalent to 

l\{2k-l)-i2k+l) a")(^"+^)(n+l)^2 ^'-'^ 
as n ^ oo. This may be established using Stirling's approximation 



V27rm i^—j 



Alternatively, there are many elementary proofs of (jl.2l) in the literature. Among 
them, [33] and [H] have recently appeared in this journal. 

Section [3] presents a proof of (|1.2I) based on the evaluation of the rational 
integral 

2 r°° dx_ 

This integral is discussed in the next section. The motivation to generalize (jl.3p 
has produced interesting links to symmetric functions from combinatorics and 
to one-loop Feynman diagrams from particle physics. The goal of this work is 
to present these connections. 



Gn z / TZTT^- (1-3) 



2 A rational integral and its trigonometric ver- 
sion 

The method of partial fractions reduces the integration of a rational function to 
an algebraic problem: the factorization of its denominator. The integral (|1.3p 
corresponds to the presence of purely imaginary poles. See [3] for a treatment 
of these ideas. 

A recurrence for Gn is obtained by writing 1 — (x^ + l) — a:^ for the numerator 
of (|1.3p and integrating by parts. The result is 

Gn+l — — Gn- (2.1) 

Zn 

Since Gi = 1 it follows that 

The choice of a new variable is one of the fundamental tools in the evaluation 
of definite integrals. The new variable, if carefully chosen, usually simplifies 
the problem or opens up unsuspected possibilities. Trigonometric changes of 
variables are considered elementary because these functions appear early in 
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the scientific training. Unfortunately, this hides the fact that this change of 
variables introduces a transcendental function with a multivalued inverse. One 
has to proceed with care. 

The change of variables x — tan 9 in the definition (|1.3p of G„ gives 

2 /""/^ 
G„+i = - / {cos9f"d9. 
Jo 

In this context, the recurrence (12. ip is obtained by writing 

(COS0)2" = (cOSf?)2"-2 + ^ii^A(cos0)2"-l 

2n — 1 dfc' 

and then integrating by parts. Yet another recurrence for G„ is obtained by a 
double-angle substitution yielding 

and a binomial expansion (observe that the odd powers of cosine integrate to 
zero). It follows that 




Thus, ()2.2[) is equivalent to proving the finite sum identity 




There are many possible ways to prove this identity. For instance, it is a perfect 
candidate for the truly 21st-century WZ method [TB] that provides automatic 
proofs; or, as pointed out by M. Hirschhorn in [TD], it is a disguised form of the 
Chu-Vandermonde identity 




(which was discovered first in 1303 by Zhu Shijie). Namely, upon employing 
Legendre's duplication formula for the gamma function 

r(i)r(2z + 1) = 22^r(z + i)r(z + i), 

the identity ()2.3|) can be rewritten as 
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This is a special case of (|2.4p . Another, particularly nice and direct, proof of 
(|2.3p . as kindly pointed out by one of the referees, is obtained from looking at 
the constant coefficient of 




Remark 2.1. The idea of double-angle reduction lies at the heart of the rational 
Landen transformations. These are polynomial maps on the coefficient of the 
integral of a rational function that preserve its value. See [H] for a survey on 
Landen transformations and open questions. 



3 A squeezing method 

In this section we employ the explicit expression for Gn given in (|2.2p to establish 
Wallis's formula This approach is also contained in Stewart's calculus text 

book [in] in the form of several guided exercises (45, 46, and 68 of Section 7.1). 
The proof is based on analyzing the integrals 

f.7r/2 



In / (sin x)" dx. 
Jo 



The formula 

follows from (|2.2[) by symmetry. Its companion integral 



'2n+l 



(2n+ 1)!! 



is of the same flavor. Here n!! = n(n — 2)(n — 4) • • • {1 or 2} denotes the double 
factorial. The ratio of these two integrals gives 

7r 

Wnl2n/l2n+l = , 

where Wn is defined by (|1.2p . The convergence of Wn to 7r/2 now follows from 
the inequalities 1 < hn/hn+i < 1 + l/(2n). The first of these inequalities is 
equivalent to hn+i < hm which holds because (sina;)^"+^ < (sinx)^". The 
second is equivalent to 

2n (sinx)2"dx < (2n+ 1) / {sinxf^+^dx, 
Jo Jo 

which follows directly from the bound /2„ < ^2n-i a-nd the recurrence {2n + 
l)J2n+i = 2n/2„_i. Alternatively, the second inequality can be proven by ob- 
serving that the function 

.7r/2 

/(s) ~ s (sinx)'' dx 
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is increasing. This may be seen from the change of variables t = sinx and a 
series expansion of the new integrand yielding 

, ^1 f2k\ 2k + 1 , , 

Remark 3.1. Comparing the series p.ip at s = with the hmit 

rTi/2 



/'(O) = lim ^ = lim / 
immediately proves 



sm a; da; = — 



^VA:/2fc+1^2' 
This value may also be obtained by letting x = ^ in the series 



E 

k=0 



2 

2fc\ x'^'^ arcsin2a; 



A: / 2fc + 1 22; 



The reader will find in [TT] a host of other interesting series that involve the 
central binomial coefficients. 



4 An example of Ramanujan and a generaliza- 
tion 

A natural generalization of Wallis's integral (jl.3p is given by 

Gn(g) = -/ li^^^-^dx, (4.1) 

where q = (gi, 52, • ■ • , 9n) with (z C. This notation will be employed through- 
out. Similarly, q" is used to denote {qf,q2, ■ • ■ , 9")- As the value of the integral 
(|4.ip is fixed under a change of sign of the parameters qk, it is assumed that 
Keqk > 0. Note that the integral Gn{q) is a symmetric function of q that 
reduces to G„ in the special case gi = • • • = 9,1 = 1 . 

The special case n = A appears as Entry 13, Chapter 13 of volume 2 of B. 
Berndt's Ramanujan's Notebooks [5], in the forn£|: 

Example 4.1. Let gi, 52, 93, and 54 be positive real numbers. Then 

2 /-^ dx 

nJo (a;2 + ql){x^ + ql)(x^ + qi){x^ + ql) ~ 

(gi + g2 + g3 + qif - (gi + gj + gj + gl) 

3glg2g3g4(gl + g2)(g2 + g3)(gl + g3)(gl + g4)(g2 + g4)(g3 + g4) ' 
minor correction from |2]. 
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Using partial fractions the following general formula for Gn{q) is obtained. 
In the next section a representation in terms of Schur functions is presented. 

Lemma 4.2. Let q — (gi, . . . , g„) be distinct and Reg^ > 0. Then 

n ^ n 

-E-n^CT- (4-2) 

k=i * j=i 

Proof. Observe first that if 6i, 62, ■ • ■ , &n are distinct then 

n ^ n 1 1 

Replacing y by a;^ and 6fc by and using the elementary integral 

2 r dx 1 



TT Jq x^ + q 

produces the desired evaluation of Gn{q)- □ 

Remark 4.3. Gn{q), as defined by (j4.ip . is a symmetric function in the q^'s 
which remains finite if two of these parameters coincide. Therefore, the factors 
Qj — Qk in the denominator of the right-hand side of (j4.2p cancel out. This 
may be checked directly by combining the summands corresponding to j and 
k. Alternatively, note that the right-hand side of (j4.2p is symmetric while the 
critical factors qj — qk in the denominator combine to form the antisymmetric 
Vandermonde determinant. Accordingly, they have to cancel. 

Example 4.4. The identities 

2 1 1 



2 poo it'^-^- -j^ 



- , . (2n + l)n!(n+ 1)!' 

j=i 

2 .oon+l ^ 1 



TT 



.^^ ^ , (2j-l)2 22"(2n+ l)(n!)2' 

2 f°" -fr 1 _ 2yl(2n- 1) 

y^^Tw or) 

may be deduced inductively from Lemma 14.21 Here, A{n, k) are the Eulerian 
numbers which count the number of permutations of n objects with exactly k 
descents. Recall that a permutation a of the n letters 1, 2, . . . , n, here written as 
cr(l) ct(2) . . . a{n), has a descent at position k if a{k) > (j{k + 1). For instance, 
A(3, 1) = 4 because there are 4 permutations of 1, 2, 3, namely 13 2, 213, 2 31, 
and 312, which have exactly one descent. 

The problem of finding an explicit formula for the numerator appearing on 
the right-hand side of (|4.2[) when put on the lowest common denominator is 
discussed in the next section. 
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5 Representation in terms of Schur functions 



The expression for Gnio) developed in this section is given in terms of Schur 
functions. The reader is referred to 4 for a motivated introduction to these 
functions in the context of alternating sign matrices and to |17j for their role 
in the representation theory of the symmetric group. Among the many equiv- 
alent definitions for Schur functions, we now recall their definition in terms of 
quotients of alternants. Using this approach, we are able to associate a Schur 
function not only to a partition but more generally to an arbitrary vector. 

Here, a vector /i — (/ii,/i2, . ■ .) means a finite sequence of real numbers. ^ 
is further called a partition (of ni) if /ii > /i2 ^ • • • and all the parts /ij are 
positive integers (summing up to ni). Write 1" for the partition with n ones, 
and denote by \{n) the partition 

\{n) = (n- l,n-2,...,l). 

Vectors and partitions may be added componentwise. In case they are of dif- 
ferent length, the shorter one is padded with zeroes. For instance, one has 
A(n + 1) = + Likewise, vectors and partitions may be multiplied by 
scalars. In particular, a ■ 1" is the partition with n a's. 

Fix n and consider q — {qi,q2, . . . ,qn)- Let fi — (/ii, /i2, . . .) be a vector of 
length at most n. The corresponding alternant is defined as the determinant 



Mi I 



I l<i.j<n 



Again, fj, is padded with zeroes if necessary. Note that the alternant a>(„) is the 
classical Vandermonde determinant 



aA(n)(9) = = n fe"*)- 

l<t.j<n 

l<z<j<n 

The Schur function Sf^ associated with the vector ^ can now be defined as 

a\(n)[q) 

If is a partition with integer entries this is a symmetric polynomial. Indeed, 
as n ranges over the partitions of m of length at most n, the Schur functions 
s^{q) form a basis for the homogeneous symmetric polynomials in q of degree 
m. 

The Schur functions include as special cases the elementary symmetric func- 
tions Ck and the complete homogeneous symmetric functions . Namely, ek{q) = 
Sifc(g) and hk{q) = S(fc)(q). 

The next result expresses the integral Gn[q) defined in (|4.ip as a quotient 
of Schur functions. 

Theorem 5.1. Let q = {qi, . . . , qn) and Req^ > 0. Then 

_ ^A(n-l)(g) _ gA(«-l)(g) , . 

sx{n+i){q) e„(q)sA(„)(q) 
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Proof. The equality e„(<7)sA(„) (q) = s^(„+i)(q) amounts to the identity 



2n-2j 




2n-27 + l 










i,3 





which follows directly by inserting the factor qi into row i of the matrix. 

From the previous definition of Schur functions, the right-hand side of (|5.ip 
becomes 

(9) (g) 
en(q) SA(n)(g) e„(q)a2A(„)(g) ' 

Observe that a2X{n){Q) — lOi" ^"'li.j = aA(n)(9^) is simply the Vandermonde 
determinant with qi replaced by qf. Next, expand the determinant a\(^n-i)+\{n) 
by the last column (which consists of I's only) to find 

aA(n-i)+A(«)(g) = e„(q) ^ aA(n-i) (g?, 92> ■ • ■ ^Qk-i^Qk+i^ • • ■ : 9^)- 

fe=i 

Therefore 

e„(g)a2A(n)(g) * ' / fj-^ 

Observe that the only terms that do not cancel in the quotient above are those 
for which i = k or j = k. The change of sign required to transform the factors 
Ik ~ 1j to qj — q\ eliminates the factor (—1)"^'''. The expression on the right- 
hand side of (|5.2p is precisely the value (14. 2p of the integral G'„(q) produced by 
partial fractions. □ 

The next example illustrates Theorem 15.11 with the principal specialization 
of the parameters q. 

Example 5.2. The special case qk — q^ produces the evaluation 

TT Jq x'^ + q-"^ q"- 1 — <7^^ 

This can be obtained inductively from Lemma 14.21 but may also be derived 
from Theorem 15.11 in combination with the evaluation (|6.ip of the principal 
specialization of Schur functions as in Theorem 7.21.2 of [18 . 

Taking the limit g 1 in (|5.3p reproduces formula (|2.2p for G„. In other 
words, (|5.3p is a g-analog [7] of (|2.2p . Similarly, 



_ °° -1 „2n 1 „2n 

= oi/4 TT t - g 1 - g 

1-^(7 li l_g2n~l 1_^2«+1 



n=l 



is a useful g-analog of Wallis's formula (11.21) which naturally appears in [5], where 
Gosper studies g-analogs of trigonometric functions (in fact, Gosper arrives at 
the above expression as a definition for tt^ while g-generalizing the reflection 
formula r(z)r(l - z) = -^). 
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The proof of Theorem 15. II extends to the following more general result. 
Lemma 5.3. 

n n , ^ 

V TT - 



where 



A = (a- 1) • A(n) - /3 • 

^ = (a-l)-A(n + l)-(/3-l)-l". 

As a consequence, one obtains the following integral evaluation which gen- 
eralizes the evaluation of Gn{q) given in Theorem 15. II 

Theorem 5.4. Let q = {qi,...,qn) and Req^ > 0. Further, let a and (3 be 

given such that a > 0, < /3 < an, and (3 is not an integer multiple of a. Then 

sm{TT(3/a) p x^-^ s^iq) 
Jo Uk^ii^ +1k) 
where A and fj, are as in Lemma \5.3[ 

Proof. Upon writing /? = 6a + /3i for 6 < n a positive integer and < /3i < a, 
the assertion follows from the partial fraction decomposition 



the integral evaluation 

f""" x^^~^dx _ 1 7r/a 



sin(7r/3i/a)' 

and Lemma 15.31 □ 



6 Schur functions in terms of SSYT 

The Schur function s\{q) associated to a partition A also admits a represen- 
tation in terms of semi-standard Young tableaux (SSYT). The reader will find 
information about this topic in [4 . Given a partition A = (Ai, A2, . . . , Am), the 
Young diagram of shape A is an array of boxes, arranged in left-justified rows, 
consisting of Ai boxes in the first row, A2 in the second row, and so on, ending 
with Am boxes in the mth row. A SSYT of shape A is a Young diagram of shape 
A in which the boxes have been filled with positive integers. These integers 
are restricted to be weakly increasing across rows (repetitions are allowed) and 
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strictly increasing down columns. From this point of view, the Schur function 
s\{q) = sx{qi, ■ ■ ■ ,qn) can be defined as 



sx 



where the sum is over all SSYT of shape A with entries from {1, 2, . . . , n}. The 
symbol is a monomial in the variables qj in which the exponent of qj is the 
number of appearances of j in T. For example, the array shown in Figure [T] is a 
tableau T for the partition (6, 4, 3, 3). The corresponding monomial is given 

by qiqlq^qlqtqlqiqs- 



1 


2 


2 


4 


5|5| 


2 


3 


4 


5 




4 


6 


6 






5 


7 


8 







Figure 1: A tableau T for the partition (6,4,3,3). 

The number Nn{\) of SSYT of shape A with entries from {1, 2, . 
be obtained by letting g — >■ 1 in the formula 



n 

l<i<j<n 



qJ-^ - 9^-1 



(see page 375 of [H]). This yields 



l<2<_;<n 



A,: - Aj + j - z 



, n} can 
(6.1) 

(6.2) 



The evaluation (12.21) of Wallis's integral (II. 3p may be recovered from here as 

^ SA(„)(1"+^) ^ Nn+i{X{n)) ^ f2n\ 
^A(„+2)(l"+') iV„+i(A(n + 2)) 22" U A 

7 A counting problem 

The fc-central binomial coefficients c{n,k), defined by the generating function 
(l-fc2x)-i/''- = Ec(n,fc)x", 



n>0 



are given by 



n-l 



c(n, k) = — - T~r (1 + km). 
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For k = 2 these coefficients reduce to tfie central binomial coefficients (^^) . The 
numbers c(n, k) are integers in general and their divisibility properties have been 
studied in [5D]. In particular, the authors establish that the fc-central binomial 
coefficients are always divisible by k and characterize their p-adic valuations. 

The next result attempts an (admittedly somewhat contrived) interpretation 
of what the numbers — c(n, —k) count. 

Corollary 7.1. Let A and fi be the partitions given by 

A = (fc-l)-A(n)-l"-\ 
M = (fc-l)-A(n + l). 

Then the integer —c{n, —k) enumerates the ratio between the total number of 
SSYT of shapes A and n with entries from {1, 2, . . . , n} times the factor A;^"~^/n. 



Proof. By Theorem 15.41 



iV„(A) sin(^/fc) r 1 , 

-G„,fc,i(l ) = — -— / , .„ dx, 



Nnifi) ' n/k Jo {x^ + l 

and the expression on the right-hand side is routine to evaluate: 

NnjX) ^ r(n-l) ^ "j^^ km -I 

N^ifi) r(n)r(i-i) }^}^ km 

This product equals —c{n,—k) divided by □ 

Remark 7.2. R. Stanley pointed out some interesting Schur function quotient 
results. See exercises 7.30 and 7.32 in 1181. 



8 An integral from Gradshteyn and Ryzhik 

It is now demonstrated how the previous results may be used to prove an integral 
evaluation found as entry 3.112 in j^. The main tool is the (dual) Jacobi-Trudi 
identity which expresses a Schur function in terms of elementary symmetric 
functions. Namely, if A is a partition such that its conjugate A' (the unique 
partition whose Young diagram, see Section [6l is obtained from the one of A by 
interchanging rows and columns) has length at most m then 

sx = \ex'^-i+j\^^^^.^^ ■ 

This identity may be found for instance in [18l Corollary 7.16.2]. 
Theorem 8.1. Let fn and gn be polynomials of the form 

gn{x) = box^^-^+bix^'''^ + ... + br,-i, 
fn{x) = aoa;" + flix""-^ + . . . + a„, 
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and assume that all roots of fn lie in the upper half-plane. Then 



9nix)dx 



-Ki Mr, 



In{x)fn{ 



where 







03 


ag . 


. 




bo 


bi 


b2 . 


• bn-i 




Oq 


02 


04 







ao 




04 










ai 


03 










ai 


03 
















an 













an 



Proof. Write fn{x) = oo 11^=1 ~ -^^y assumption, Reg^- > 0. Further, 



fnix)fni'x)^i-iralllix^ + q^). 

Let q = {qi,q2, ■ . ■ ,qn)- It follows from Theorem 15.41 that 

x^^dx _ (-l)"+'^7rs^(„_i)_2/3.i"-i(g) _ (-l)"7r sx'iq) 



— CO fn{x)fn{ x) 



S\(n+1)-2P-1" (9) 



ao SA(„+i)(g) 



where A = X{n — 1) + 2 • l'^. The latter equality is obtained by writing the 
quotient of Schur functions as a quotient of alternants, multiplying the fcth row 

2B 

of each matrix by , and reordering the columns of the determinant in the 
numerator. The right-hand side now is a quotient of Schur functions to which 
the Jacobi-Trudi identity may be applied. In the denominator, this gives 

S\{n+l){q) = \en+l-2k+]{q)\i<k,j<n = I ^2*:-^ (?) I i<fc j <„ • 

Note that em{q) = whenever m < or m > ?i. Further, ek(q) = i'^ak/ao. 
Hence, sa()i+i)('7) = i"'"^^-'^^A„/aQ . The term sx'{q) is dealt with analogously. 
The claim follows by expanding the determinant Af„ with respect to the first 
row. □ 



9 A sum related to Feynman diagrams 

Particle scattering in quantum field theory is usually described in terms of Feyn- 
man diagrams. A Feynman diagram is a graphical representation of a particular 
term arising in the expansion of the relevant quantum-mechanical scattering am- 
plitude as a power series in the coupling constants that parametrize the strengths 
of the interactions. 

From the mathematical point of view, a Feynman diagram is a graph to 
which a certain function is associated. If the graph has circuits {loops, in the 
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physics terminology) then this function is defined in terms of a number of in- 
tegrals over the 4-dimensional momentum space {ko,k), where kg is the energy 
integration variable and fc is a 3-dimensional momentum variable. 

Feynman diagrams also appear in calculations of the thermodynamic prop- 
erties of a system described by quantum fields. In this context, the integral over 
the energy component of a Feynman loop diagram is replaced by a summation 
over discrete energy values. These Matsubara sums were introduced in |14| . A 
general method to compute these sums in terms of an associated integral was 
presented in [S]. 

These techniques, applied to the expression (|4.2p for the integral G'„(q), give 
the value of the sum associated with the one-loop Feynman diagram consisting 
of n vertices and vanishing external momenta, Ni = 0, as depicted in Figure [2] 




2 



Figure 2: The one-loop Feynman diagram with n vertices and vanishing external 
momenta, m is the summation variable associated to each of the internal lines. 



The Matsubara sum associated to the diagram in Figure [2] is 



oo n 



J. X _|_ 



m— — OO k—1 



where the variables qk are related to the kinematic energies carried by the (vir- 
tual) particles in the Feynman diagram. This sum was denoted by Sg in [5]; 
the notation has been changed here to avoid confusion. 

Example 9.1. The first few Matsubara sums are 
91 

. , , X q2Di - qiD2 
M2{qi,q2j = T^—r-n 2T' 



^3(91,92,93) = TT 



9192(95 -91)' 

9293(9! - 9i)-Di + 9391(9;? - ql)D2 + 9192(9? - 9i)£'3 



919293(9! - 9i)(9l - 9?)(9? - 9I) 
with Dj = coth(7rgj). 
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Theorem 9.2. The Matsubara sum Mn{q) is given by 



n-T r \ coth(7rgfe) -A- 1 



k=l 



Proof 1. This follows from the partial fraction expansion 
" 1 " 1 1 

rr t = y \ tt ^ 



k=l 



which is a special case of (j4.3p . switching the order of summation, and employing 
the classical 

TT coth(7rz) 

z2 + w? 



00 



Proof 2. The method developed in 5 shows that 



M„(q-) = TT 



1 + ^ nb(gm)(l - Rra) 



Gniq), 



□ 



(9.2) 



where Gn{q) is the integral defined in (|4.ip . 



nb((7) 



1 



1 



- (cothTrg - 1) , 



g27rq _ 2 2 

and i?m is the reflection operator defined by 

Rmfiqi, . . . , g™, . . .) = f{qi, -qni, ■ ■ •)■ 

To use (|9.2I) combined with the evaluation (I4.2p of G„ (q) it is required to com- 
pute the action of each 1 — R,„ on the summands of (|4.2p . Namely, 



n 

(i-i?™)-n— 



2(5, 



km 



n ^ 

n „2 _ „2 ' 



where Sh„, is the Kronecker delta. Therefore 



nb(gm)(l - i?m)G„(q) 



2nb(gn 



' n „2 „ 2 



and the result follows from 2nb(q) = coth{TTq) — 1. 



□ 
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Finally, an expansion of Mn{q) in terms of symmetric functions is given. 
Starting with the classical expansion 

where C(s) denotes the Riemann zeta function, it follows that 

71 OO Tl 

M..(.) ^ E 4 n ^ - 2 E (-i)"<(^".) i: n 

k=l " j^k J " m=l fe=l j^k J ^ 

Using the identity {hj being the complete homogeneous symmetric function) 

n 

fe=l ]^k ■' 

which follows from Lemma [573l for see page 450, Exercise 7.4 of [18]), this proves: 
Corollary 9.3. The Matsubara sum Mn{q), defined in Ii9.1\) . is given by 

Mn{q) = + 2 (-l)"C(2m + 2n)h^{q^). 

^"(.9 I m=0 

10 Conclusions 

The evaluation of definite integrals has the charming quality of taking the reader 
for a tour of many parts of mathematics. An innocent-looking generalization of 
one of the oldest formulas in analysis has been shown to connect the work of 
the four authors in the title. 
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